Coexistence of Enhanced Superconductivity and Antiferromagnetism: Possible Correlated Phase 

Transitions in Trilayer High-T c Cuprates 
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Based on a hybrid inteiiayer coupling mechanism, we study the coexistence of superconductivity (SC) and 
antiferromagnetism (AFM) in trilayer cuprates. By introducing an interlayer magnetic scattering term, we solve 
the multilayer t-J model with Josephson coupling under the framework of Gutzwiller projection. We show that 
both the SC and AFM orders in the multilayered system are enhanced and the range of AFM order is extended. 
The layer configuration of rf-wave pairing gap and AFM order further plays an essential role in determining the 
interlayer magnetic and superconducting coupling phase diagram of such multilayered systems. Abrupt phase 
transitions between correlated states carrying distinct configurational symmetries are unveiled by tuning the 
doping level and/or the tunneling strengths. 

PACS numbers: 74.78.Fk, 74.72.Gh, 74.25.Dw 



Motivation. — Strong electron correlations and unconven- 
tional order parameters fuse together to engender a plethora 
of fascinating phenomena in multilayered high-T c cuprate su- 
perconductors (HTSC). Among them is the discovery of the 
coexistence of HTSC phase with the AFM metallic phase in 
site-selective nuclear magnetic resonance (NMR) studies of 
multilayer cuprates liHSD- Experimentally, the presence of 
outer CuC<2 layers provides the effective screening of the ran- 
dom potential influence from the charge reservoirs in the in- 
ner CuC<2 planes. This intrinsic feature greatly facilitates the 
realization of the coexistence of HTSC with AFM in multilay- 
ered cuprates and makes their mutual interplay a central issue 
as compared to single layered systems. Furthermore, a strik- 
ing common aspect shared by the phase diagrams of a large 
class of SC materials including cuprates, iron-pnictides and 
several heavy-fermion superconductors is the close proximity 
of SC phase with the AFM or spin-density- wave phase @]. 
Therefore, such an intriguing mixture of HTSC and AFM in 
multilayered copper oxides sets up a new platform for investi- 
gating the relationship between these two orders. It also opens 
the key avenue toward uncovering the mechanism underlying 
HTSC in cuprates and related materials. 

Angle-resolved photoemission spectroscopy (ARPES) 
measurements on the high-quality optimally doped triple- 
layer cuprate superconductor Bi2Sr2Ca2Cu30io+<5 (Bi2223) 
further reveal a layer variation of both doping density and d- 
wave SC gap through successfully observing the electronic 
structure and multilayer band splitting @]. In particular, two 
symmetric outer CuC>2 planes (OPs) are overdoped with gaps 
which are larger than those for optimally doped single layer 
cuprates while the inner CuC<2 plane (IP) is underdoped with 
an even larger gap. Generally, an interlayer Cooper pair tun- 
neling has always been put forward as an approach to increase 
the SC critical temperature of multilayer cuprates @, Q] and 
was extensively elaborated to explain why T c takes maximum 



at ?i = 3 Nevertheless, there is still a significant ques- 
tion concerning the role of Cooper pairs tunneling in the AFM 
background when the uniformly coexistent AFM-SC state was 
discovered in multilayered copper oxides. 

Theoretically, although it is known that two-dimensional 
t-J model allows the emergence of AFM order and its coexis- 
tence with HTSC in the underdoped region jg], the underlying 
mechanism of enlarging the coexisting range between AFM 
and HTSC phases in multilayer cuprates is still under con- 
troversial debate |I3-i5]. To this end, a proper understand- 
ing of the exhibited broadening of the AFM-SC mixed regime 
remains much needed. Moreover, both SC and AFM orders 
carry nontrivial phase factors in the geometrically nonequiv- 
alent IP and OPs, and the detailed structure of the phase di- 
agram of the multilayered cuprates constitutes an important 
problem for study. However, all these issues have not yet been 
fully investigated so far in the literature. 

In this Letter, we study the microscopic coexistence and 
interplay between HTSC and AFM orders within a proxim- 
ity scenario. An exotic intramultilayer AFM scattering term 
is proposed to elucidate the appreciable enlargement of the 
range for the mixed phase. Besides, once isolated Q1O2 
planes coupled together by interlayer tunneling processes, 
originally degenerate states will be broken down to certain 
peculiar configurations, which allows us to identify the sig- 
nificance of such layer arrangements of SC and AFM phases 
in settling the ground state of multilayered cuprates. 

Multilayered t-J model and Gutzwiller projection. — In the 
present paper we shall focus mainly on the system of trilayer 
HTSC copper oxides and treat it as a prototype to probe the 
SC and AFM properties of multilayer cuprates. For simplic- 
ity, the possibilities of interlayer electron pairing and direct 
interactions between two OPs are ignored, and we only con- 
sider intralayer Cooper pairing and adjacent interlayer cou- 
plings within one unit cell of a trilayer cuprate. t-J model 
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is adopted to describe the microscopic electronic structure of 
each layer. The effective free energy function describing such 
trilayer systems can be written as 

a</3,<r^ ^ ' 

-2 3s z Q Jm2-2^r^-4 7 A Q A /3 +C/ J _ ( S Q d>- MQ n Q }(l) 

per site, where three CUO2 planes (layer-index a, (3 — 1, 2, 3) 
share the same parameters of t-J model and the same inter- 
layer Coulomb repulsion U± iflrjll . which could redistribute 
the doped holes over individual layers. A phenomenological 
pair tunneling Hamiltonian originating from a Josephson-like 
pair tunneling process is included in T lfl7ll . Two nearest- 
neighbor electrons with opposite spins are considered to con- 
dense into a Cooper pair tunneling in the real space. In 
Equation (Q~|), we define such singlet pairing amplitudes as 
A Q = (cf^c^ — c[°^c^)o- The notation (• • -)o indicates 
the expectation value in the unprojected state. The other two 
order parameters decoupled from the spin-spin superexchange 
interaction are defined as Xa = ( c if c jf 



C i,i C j,l>° aIld 

m a = e'"<"' (Sf Q )o- represents a spin-oriented magnetic 
scattering between adjacent layers deduced from the channel 
of single electron tunneling, whose form is as follows [18]: 
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FIG. 1: (color online). The staggered AFM moments m a and the 
projected SC order parameters Asc.a in the IP and OPs for the 
trilayer system as functions of individual doping levels 8 a , with 
K—0.23t, 7=0. 13t, U±=1.6t. Corresponding orderings m and Asc 
for the single-layer t-J model were shown for comparison [panel a]. 
The carrier distributions in IP and OPs are plotted as functions of 
average doping S [panel b]. The calculated AFM moment versus 
doping in each layer with the chosen parameters was compared with 
cases of k=7=(7L=0 (m, black curves) and /t=0, 7=0. 13t, U±=1.6t 
(m' a , light red and blue curves) for IP [panel c] and OPs [panel d]. 
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The summation of k runs over the reduced Brillouin zone: 
— 7r < k x + k y ^ 7T, the nesting vector Q = (tt^jt), and 
^k = cos k x — cos k y is a <i-wave symmetrical factor 11911 . Here 
k and 7 are strengths of interlayer AFM scattering and Cooper 
pair tunneling. By deploying extended Gutzwiller approxi- 



mation 11201 12111 to relax the constraint of no doubly occupied 
sites, we substitute the projected order parameters with the 
above-defined mean fields in the unprojected Hilbert space 
multiplied by various statistical counting factors, which re- 
sults in the free energy function J 7 . The detailed analytic 
forms of these Gutzwiller factors can be found in Ref. 1 2111 
and a simplified version for the renormalized factor g"^_ is 
effectively expressed in terms of the doping levels and AFM 
moments of IP and OPs OH, 
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This novel interlayer scattering between different magnetic 
Brillouin zones mimics a hidden AFM correlation in the intra- 
multilayer magnetic structure, and it produces a desired con- 
sequence of enhancing AFM moments with the HTSC intact. 
Apparently, such kinds of interlayer magnetic couplings are 
solely arising from the strong electron correlation effects in 



underdoped cuprates and once the AFM moments vanish, this 
channel will be switched off owing to the shrinking of j_. 
Throughout the whole numerical calculation, we set t = 300 
meV as the energy unit and J — 0.3t, which are typical pa- 
rameters in the t-J model. 

Enlarged coexistent region of AFM and HTSC. — Following 
standard variational procedure, we minimize the free en- 
ergy T to obtain varied order parameters from a set of self- 
consistency equations. The stable local energy minima can be 
found within restricted parameter ranges by iteration method. 
In the presence of commensurate AFM order, the Brillouin 
zone is fold into two parts, which enables us to construct intra- 
and interlayer magnetic scatterings with nesting vector Q. 
The fully entangled Gutzwiller renormalization factors further 
generate strongly nonlinear effects on the self-consistently de- 
termined mean-field order parameters, and the explicit depen- 
dence of g" j_ on the spin indices and AFM moments from 
nonequivalent IP and OPs plays another essential role to lift 
the degeneracy of intramultilayer magnetic structures. 

In Fig. [it, we plot the staggered magnetizations m a and 
the projected SC order parameters Asc, a = .9aA q in both 
IP (red curves) and OP (blue curves) as functions of individ- 
ual dopings 8 a . The solid and dash lines represent AFM mo- 
ments and projected SC gaps, respectively. For comparison, 
we have used the black solid and dash lines to denote the re- 
sults of decoupled single-layer t-J model. A critical doping 
S c can be resolved around 0.1 in the black solid line, below 
which AFM order coexists uniformly with SC order 12111 . As 
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the values of parameters increasing to k = 0.234, 7 = 0.134 
and U±=1.6t, the ground-state layer configuration of the tri- 
layer system in the mixed regime favors an intramultilayer 
AFM arrangement with coherent superconductivity. The stag- 
gered magnetic phases m a also become much stiffer against 
dopings for both inner and outer layers. The resulting quan- 
tum critical doping points <5 C , Q beyond which m a diminish 
steeply are 0.14 and 0.17 for the IP and OPs, respectively. 
Our result qualitatively agrees with the NMR observation of 
trilayer cuprates Ba2Ca2Cu30e(F,0)2 0]. This increase of 
8 c . a to higher magnitudes arises from interlayer AFM scat- 
terings and the high entanglement of order parameters in the 
limit of strong electron correlation. The interlayer quantum 
tunneling of Cooper pairs further immensely promotes the SC 
gaps throughout the whole system, and the maximum value 
of the projected SC order parameter in the IP is nearly twice 
larger than that of the single-layer cuprate, which is consis- 
tent with the enhanced superconducting gap observed by the 
ARPES experiment in Bi2223 @]. Moreover, around S c , Q , the 
SC orders for both IP and OP always reach their maxima. 

Another subtle phenomenon in realistic multilayered mate- 
rials is the nonhomogeneous doping distribution among inner 
and outer cuprate planes jll-blpL 12211 . We describe such an es- 



sential effect by U± lfl6ll . Considering the doped holes have a 
gradually diffusive process from outer planes to the inner, two 
different values have further been assigned to {fi a } to exem- 
plify a possible layer gradient of chemical potential. The self- 
consistently calculated charge distribution is shown in Fig.rjJ) 
with the same values of parameters as in panel a. We find that 
throughout the average doping axis, the doping level of OP is 
always greater than that for the IP and their imbalanced value 
is about 0.05 when the averaged doping S equals 0.19. This is 
consistent with the experimental observations 12211 . Further- 
more, the concentrations of doped holes in the inner and outer 
planes are both monotonically increasing functions of S, ex- 
cept that there exist kinks at the critical average doping 0.15, 
beyond which the system is purely superconducting. Finally, 
the coexistence of AFM and HTSC might slightly shrink the 
resulting charge imbalance between IP and OP. 

Due to the competition between emergent AFM and en- 
hanced HTSC, the values of AFM moments of the trilayer 
system would be greatly suppressed in the region close to 
S c , if the scattering strength k was reduced to zero. How- 
ever, once K becomes comparable to 7, this suppression can 
be reversely changed, especially in the regime with large dop- 
ing concentrations, as demonstrated in Fig. [TJ; and d, where 
we compare the results of staggered magnetization for cases 
with sets of parameters: k = 0.234, 7=0.134, U± = 1.6t and 
K=0, 7=0. 13t, U± = 1.6t, as well as k=j = U±=0 in both 
IP and OP. One of the unusual features of the designed in- 
terlayer AFM coupling is its compatibility with the interlayer 
Cooper pair tunneling, which means we can conversely com- 
pensate the suppression of AFM moments via magnetic scat- 
tering without reducing the enhanced superconductivity from 
the channel of pair tunneling. This unique combination and 
interplay between these two channels underpins the realiza- 



tion of the extended coexistence of AFM and HTSC in our 
trilayer t-J model and the achieved charge imbalance further 
facilitates this phenomenon. 

Predicted phase diagrams of trilayer cuprates. — We would 
like now to analyse the relative layer configurations of coexist- 
ing AFM and SC order parameters in the presence of both in- 
terlayer magnetic scattering and Cooper pair tunneling. Since 
these two orderings have the Ising symmetry, the interlayer 
couplings will drive four distinct stacking patterns between 
IP and OP in one unit cell when only thermodynamic stable 
phases are considered. A candidate state here is specified by 
its phase arrangements of both AFM and SC orders among 
three layers. We distinguish the layer arrangements of AFM 
moments into two categories: the ferro-type structure, namely 
[sgn. (™"p =l], which indicates that the interlayer configura- 
tion of AFM moments is ferromagnetically aligned in one unit 
cell, and the antiferro-type structure, namely [sgn(^^-)=— 1], 
denoting the AFM arrangement. In the same way, we have in- 
phase [sgn(|^) = l] and out-of-phase [sgn(^) = -l] SC 
states. Therefore, totally there exist four possible degener- 
ate ground states with distinct symmetries in layer configu- 
rations of order parameters when three layers are decoupled 
(/t=7=0). However, interlayer tunnelings, as weak pertur- 
bations, will spoil originally degenerate equilibrium state to 
some special symmetry broken configuration. Possible abrupt 
phase transitions between these correlated states will be un- 
cloaked through tuning the strengths of n and 7 and/or the 
averaged doping S. In the following, we will use (±,±) to 
denote these four states individually and the first and second 
"±"s represent the intramultilayer arrangements of AFM and 
SC orderings, respectively. The first "+" sign on the left refers 
to the ferro-type structure for the AFM order and the second 
"+" sign on the right refers to the inphase state for the SC or- 
der, while the "— " signs are adopted for the opposite magnetic 
and superconducting configurations. 

Minimization of the free energy T in the order parametric 
space leads to the genuine ground state with specified config- 
urational symmetries. The internal competition, entanglement 
and renormalization of orderings further combine together to 
manifest the essential interplay between AFM and HTSC in 
multilayer materials. Physically, we know that sufficiently 
large interlayer magnetic scattering will drive the triple-layer 
system to exhibit the antiferro-type magnetic structure, while 
the strong Cooper pair tunneling will give rise to the inphase 
SC state. Nevertheless, the detailed evolution, when sweep- 
ing through the relevant parametric space (k, 7, 5), is still not 
clear for multilayered t-J model. In view of the qualitatively 
consistent results derived above from the trilayer proximity 
model, in the following, we will address this crucial problem 
based on analysing the same free energy function T . 

We first consider the weak Josephson-coupling limit, from 
which we can extract some generic features about the tran- 
sitions between these correlated ground-state configurations. 
Figure [2^ is a phase diagram for the AFM and HTSC coex- 
isting region plotted in terms of k and S with 7 = 0.044 and 
U± fixed to be 1.64. Four possible layer configurations are all 
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FIG. 2: (color online). The ground-state phase diagram for the triple- 
layer system in the weak Josephson-coupling limit with 7 = 0.04t 
[panel a], as well as the phase diagram in the strong limit of Joseph- 
son coupling with 7 = 0. If [panel b]. The inset of panel b depicts 
the free energy difference between states (+, +) and (— , +) when 
tuning k at 5=0.065. See text for the definitions of states (±, ±). 



present in different parts of the k-6 plane. In the almost un- 
doped regime, the antiferro-inphase state (— , +) always pos- 
sesses the lowest free energy than other arrangements along 
the k axis. When gradually increasing 5 to the intermedi- 
ate underdoped regime, the antiphase SC state emerges as 
the ground state with a variety of magnetic structures: one 
first enters a region in which the symmetry (— , — ) domi- 
nates, and then with increased k, state (+, — ) supersedes 
( — , — ) to be energy favored. As further lifting k above a 
critical value, the system will undergo an abrupt phase tran- 
sition accompanying with both SC and AFM phase-changes: 
(+, — ) (— , +). These entangled simultaneous phase-shifts 
offer a unique manifestation of the interplay and correlation 
between the symmetries of AFM and SC phase arrangements, 
and the sufficient renormalization of chemical potential, the 
resulting nonhomogeneous charge distribution and the energy 
competition between AFM, SC, and kinetic orderings are the 
major underlying driving force of this transition. Experimen- 
tal route to probe the signature of such a phase transition might 
include using phase-sensitive facilities 12311 . As S further ap- 
proaching 5 C , the intramultilayer SC phases in one unit cell 
will incline to form a coherent arrangement with a phase tran- 
sition from the ferro-type magnetic structure to the antiferro- 
type when continuously promoting the strength of k. Beyond 
5 C , the AFM ordering of the system is mainly sustained by 
the channel of interlayer magnetic scattering, thus the only 
allowed magnetic structure is of the antiferro-type. 

In the strong limit of interlayer Josephson coupling, the out- 
of-phase SC state will become unstable and eventually fade 
out for all the doping concentrations we consider. The inphase 
SC state prevails to be the ground state with various magnetic 
structures when 7 = O.lt. We find that if the averaged dop- 



ing is within intermediate underdoped regime 0.06<<5<0.09, 
the ferro-inphase state (+, +) will have a lower free energy 
than the antiferro-inphase state (— , +) when k is varying be- 
tween Ot to 0.2t. The resulting phase diagram is presented 
in Figure [2J5, where we calculate a dome-like phase bound- 
ary separating these two phases. To gain a concrete under- 
standing of this magnetic phase transition, we plot the free 
energy difference between states (+, +) and (— , +) as a func- 
tion of k in the inset of Figure ^p. The average doping is 
chosen to be 0.065 and the energy difference is defined as 
AJ 7 = ^(+.+) — J r (_.+). The boundary between these two 
phases can be obtained by nontrivial values of k and 5 sat- 
isfying AJ 7 = 0. In the other region of the phase diagram, 
the ground state of the trilayer system should have the (— , +) 
symmetry and exhibit the antiferro-type magnetic structure 
with enhanced superconductivity. 

While our minimal model (Q~|i is simple, implications of the 
AFM-SC mixed phase and the mapped phase diagrams can be 
readily generalized to more complex layer materials. Specifi- 
cally, if we increase n from 3 to 4 and to 5, the system might 
turn out to decay from a fully inphase SC state to the partially 
out-of-phase SC state with narrower pairing gaps, thus giving 
rise to a lower T c . In addition, with more layers stacking to- 
gether in one building block of multilayered cuprate, the layer 
gradient of chemical potential will be amplified in the influ- 
ence of interlayer magnetic scattering, and we need to dope 
more holes into the system to completely suppress the AFM 
order in the innermost CUO2 plane, which reversely indicates 
the microscopic coexistence of AFM with HTSC can survive 
to a larger average doping. 

Conclusion. — In summary, we have performed a system- 
atic investigation of the interlayer coupling effects on trilayer 
cuprates, which brings into full focus the issues concern- 
ing novel intramultilayer AFM scatterings in producing the 
observed broadening of AFM and HTSC coexisting region 
and abrupt correlated phase transitions between four different 
AFM-SC states for the trilayer system. Our predicted phase 
diagrams might serve as the first attempt to clarify symmetry- 
sensitive structures in layer materials and also to highlight the 
profound consequences of interlayer phase arrangements in 
triggering the unconventional interplay between HTSC and 
AFM in multilayered high-T^ superconductors. 
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DERIVATIONS OF THE INTRAMULTILAYER ANTIFERROMAGNETIC SCATTERING AMPLITUDE V af) AND 
THE GUTZWILLER RENORMALIZATION FACTOR g"f ± 

In this supplemental material, we will first provide the motivation for introducing the interlayer antiferromagnetic scattering 
and then give out detailed deductions for the explicit forms of and g" ± - 

To avoid possible confusions, we need to define and fix two conventions when dealing with antiferromagnetism in the extended 
Gutzwiller approximation of a single-layer t-J model in a square lattice containing N sites: 

Q_; %-J = -* E (ic j ,+H.c.) + J^S I -S j -/ i ^fi i , (1) 

(i,j),<r a,i 

,i_S ■ where c' and c represent the electron creation and annihilation operators in the projected Hilbert space. In the presence of anti- 
J]3 | ferromagnetism, we shall divide the square lattice into two sublattices: sites A with net spin up and sites B with net spin down. 
Q Under the Gutzwiller projection, we substitute each operator c\ a (ci a ) in the hopping terms with a corresponding unprojected 
operator c\ a (ci a ) multiplied by a classical Gutzwiller factor g\ a . For clarity, we demonstrate the first convention we use in 
defining a simplified version of g\ a as follows: 
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where 8 is the doping concentration and m = e ?< ^ Ri (Sf }o is the staggered magnetic moment in the unprojected state (2), When 
Fourier transforming the Hamiltonian from one space to the other, we adopt the following correspondence between operators 
CAka, CBka and c k ,<T, c k +Q,o- & 

CKa = — 71 — ' Ck+Q "~ — 71 — ' () 

Armed with these conventions, we now illustrate straightforwardly that for a multilayered t-J model, the proposed amplitude 
of the spin-oriented magnetic scattering can be regarded as an effective order parameter of the intramultilayer magnetic 
structure. Temporarily, we neglect the d-wave factor ?y£ and redefine T% as 

t° a - — V/r (Q)t c (,3) -c (Q)t c (W +Hc\ (4) 

k 

where the summation of k is summed over the magnetic Brillouin zone. It is easy to show that if the AFM moments of IP (a) 
and OP ((3) are aligned ferromagnetically in one unit cell, namely sgn(— ) = 1, the value of T% is close to zero: 
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On the contrary, when the staggered AFM moments are inclined to be antiferromagnetically arranged in a unit cell, namely 
sgn( ^ ) = — 1, the value of T^p is approaching a finite number: 
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W«2V«+H.c. - E^MJ+Kc. 1)^0, (6) 




owing to the emergence of antiferromagnetic ordering and the Pauli exclusion principle. After averaging over the d-wave 
symmetrical factor ?/ k , we find that under the same set of parameters, numerically the value of calculated in the antiferro- 
type magnetic structure is one order greater than the value of evaluated in the ferro-type magnetic structure. Motivated by 
the above observation, we dub T^p and the related process a sort of intramultilayer antiferromagnetic scattering (coupling). 

With the identification of Y^a as the amplitude of an interlayer magnetic scattering, we are going to construct an across-the- 
interface interaction which has an analogous analytic form like rjTg. The methodology we deploy here is instead of direct 
model building, we choose to filter out such kinds of magnetic interactions from the various channels of interlayer single 
electron tunneling. Within the filtering process, we will find the explicit form of the Gutzwiller renormalized factor ± . 
As an illustration, we first consider the onsite single electron tunneling between adjacent IP (a) and OP (j3) in one unit cell of a 
trilayer cuprate which is exhibiting the antiferro-type magnetic structure: 
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Under the extended Gutzwiller approximation, the renormalized Hamiltonian of W± is 

__ t V / *(") J Q )t J/3) , *(«)*(/?) JWJ") t(a) t(/3) (a)f (0) t(a) t(fi) (/3)f (a) 

n -L — [9 At 9 At C At C At 9 At 9 At c At C At 9 At 9 At c Al c Al 9 A\ 9AI c Al c Al 



*(") *(/3)J«)t J/3) , t(a) *(/?) (/3)t (a) t(a) t(/9) (a)f 03) t(a) t(P)(f))t (a)\ m 
i/Bt »St B T Bf "St fBt Bt C Bt ' 9BI 9BI c Bl c Bl ' 9b\ 9bI c Bl c Bl ) ' W 



where we have adopted the simplified version of Gutzwiller projection factor g\ a defined in Equation (0. Since we have assumed 
working in an an 
of Equation (f2]): 



working in an antiferro-type magnetic structure, we shall clearly distinguish g^. and g 1 ^ to be consistent with the convention 
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Here we have classified the sites with net spin orientation up in the IP (a) into the sublattice A and those with net spin orientation 
down in the IP into the sublattice B. By using the Fourier transformations: 
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we express H± in the momentum space as 
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where a <i-wave symmetric factor is introduced ij_(k) = ^(cos/s^. — cosk y ) 2 j^]. However, once we translate the above 
expression into the k and k+Q formalism, the desired new terms will emerge due to the special form of g\ a . According to the 
following relations: 
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(13) 

finally, we find that the renormalized Hamiltonian H± of the single electron tunneling can be cast in the form 
H± = — tj_(k) {gm9m~9A^9Mj f c k.t^ c k+Q.t~ c k+Q,t c kA+H- c ' 

k 
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where t±(k.) is modified by a factor of one half and g\^"\ g\^ are defined in Equations (0 and 1101 . respectively. Several 
comments should be noticed: 

(i) Owing to the discrepancy of g* CT in different sublattices, novel terms whose forms are identical to that of the amplitude of 
intramultilayer antiferromagnetic scattering V^ g have been emerging in the first two rows, which explicitly indicates the 
existence of possible interlayer antiferromagnetic correlations in multiplane cuprates; 

(ii) From Equation ( fT4b . we can directly read out the detailed expression for the effective Gutzwiller renormalization factor 

a/3 _ t(a) t(/3) t(a) t{B) 
9o,± — 9 Bo 9 Bo 9Ao 9Ao ' 

(iii) Because the renormalized factor g\ a is depending on the spin orientation, originally equally-weighted channels of spin up 
and spin down are now showing slightly different behaviors. In other words, based on the extended Gutzwiller approxima- 
tion, two novel channels become emergent and the resulting four channels are renormalized by g\ a with lightly different 
weights; 

(iv) We have checked that the terms in the last two rows in Equation (TBI) correspond to a Cooper-pair breaker, and have no 
direct relations with the magnetic properties of the trilayer system; 

(v) Since for our free energy function T, we allow both ferro- and antiferro-type magnetic structures in V a aB , it shall be 
justified that there exist other identical channels of the proposed intramultilayer antiferromagnetic scattering when the 
AFM moments in one unit cell is ferromagnetically aligned. A quick calculation shows that the nearest-neighbor interlayer 
tunneling of single electron just provides such equivalent channels of intramultilayer AFM coupling if we start with 
sgn(^j) = 1, therefore we can effectively use the strength k to denote all these interactions from the onsite to the nearest- 
neighbor single electron tunnelings between adjacent copper oxide layers and relax the magnetic structure to be either 
ferro- or antiferro-types, namely sgn(^-) = ±1, in T^. This finalizes the form of the trilayer proximity model we 
studied; 

(vi) Normally, the ferro-type magnetic structure is with suppressed AFM moments, whilst the coherent SC state carrying an 
inphase arrangement could greatly facilitate the <i-wave pairing amongst all three layers in collaboration with an effective 
Cooper pair tunneling. 
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We shall emphasize that all the above consequences are solely stemming from the strong electron correlation effects in under- 
doped high-r f cuprate superconductors. 
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